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Abstract 

Non-DD decay ^(3770) — > p7r is searched for using a data sample of (17.3 ± 0.5) p6 _1 taken 
at the center-of-mass energy of 3.773 GeV by the BESII detector at the BEPC. No pir signal is 
observed, and the upper limit of the cross section is measured to be a(e + e~ — > pir) < 6.0 pb at 90% 
C. L. Considering the interference between the continuum amplitude and the ^(3770) resonance 
amplitude, the branching fraction of ^(3770) decays to pir is determined to be £>(V'(3770) — > pir) € 
(6.0 x 10~ 6 , 2.4 x 10~ 3 ) at 90% C. L. This is in agreement with the prediction of the S- and 
.D-wave mixing scheme of the charmonium states for solving the "pir puzzle" between J/ijj and 
if)(2S) decays. 

PACS numbers: 13.25.Gv, 12.38.Qk, 14.40.Gx 
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I. INTRODUCTION 



A. "p7r puzzle" and -0(3770) — »■ pit 

Perturbative QCD (pQCD) predicts that the decays of J ftp and ip(2S) (shortened as ip' 
below) into light hadrons are dominated by the annihilation of cc into three gluons, with 
widths proportional to the square of the wave function at the origin [1]. This yields the 
pQCD "12% rule", that is 

Q h = p!=±- = „ 12%. (l) 

The violation of the above rule was first observed in pn and K* + K~ +c.c. modes by Mark 
II [2] , known as the pn puzzle. Since then BES and CLEO-c have measured many two-body 
decay modes of ■?/>', among which some obey the 12% rule while others violate it [3, 4]. There 
have been many theoretical efforts trying to solve the puzzle [5] . A recent one is the S- and 
.D-wave charmonia mixing model proposed by Rosner [6]. In this scheme, the mixing of the 
ip(2 3 Si) and ip(l 3 Di) states is in such a way that there is almost a complete cancellation 
of the decay amplitude of ip' — > pn, which instead shows up as an enhanced decay mode of 
■0(3770) (shortened as ip" below). A study shows that in e + e" experiments with B(ip" — > pn) 
predicted by the S- and D-wave mixing in Ref. [6], the destructive interference between the 
three-gluon decay amplitude of the ip" resonance and the continuum one-photon amplitude 
leads to a very small pn cross section at the ip" peak [7], which is in agreement with the 
unpublished upper limit of the pn cross section at the ip" peak by Mark III [8] . Using a larger 
data sample to further study this channel will shed light on the understanding of the u pn 
puzzle" and the interference pattern between the resonance and the continuum amplitudes. 



B. Data samples and detector 

The data used for this analysis are taken with the BESII detector at the BEPC storage 
ring at the center-of-mass energy of 3.773 GeV. The integral luminosity of the data sample 
is (17.3 ± 0.5) pb^ 1 as measured using large angle Bhabha events. To study the continuum 
process, BESII also collected (6.42 ± 0.24) pb' 1 data at y/s = 3.65 GeV [9]. We analyze 
these two data samples simultaneously to get the numbers of events at the ip" peak and at 
the continuum. 
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BESII is a conventional solenoidal magnet detector that is described in detail in 
Refs. [10, 11]. A 12-layer vertex chamber (VC) surrounding the beam pipe provides trigger 
information. A 40-layer main drift chamber (MDC), located radially outside the VC, pro- 
vides trajectory and energy loss (dE/dx) information for tracks over 85% of the total solid 
angle. The momentum resolution is a p /p = 0.017^1 + p 2 (p in GeV/c), and the dE/dx 
resolution for hadron tracks is ~ 8%. An array of 48 scintillation counters surrounding the 
MDC measures the time-of- flight (TOF) of tracks with a resolution of ~ 200 ps for hadrons. 
Radially outside the TOF system is a 12 radiation length, lead-gas barrel shower counter 
(BSC). This measures the energies of electrons and photons over ~ 80% of the total solid 
angle with an energy resolution of ge/E — 22%/\fE (E in GeV). Outside of the solenoidal 
coil, which provides a 0.4 Tesla magnetic field over the tracking volume, is an iron flux return 
that is instrumented with three double layers of counters that identify muons of momentum 
greater than 0.5 GeV/c. 

C. Monte Carlo Simulation 

A Monte Carlo simulation is used for the determinations of the mass resolutions and 
the detection efficiencies. This program (SIMBES), which is Geant3 based, simulates the 
detector response, including the interactions of secondary particles with the detector mate- 
rial. Reasonable agreement between data and Monte Carlo simulation has been observed in 
various channels tested [12], including e + e~ — > ( y r y)e + e~, e + e~ — > {j)p + p~ , J/ip — > pp and 
ip' -> J/^TT+vr", J/ip -> £ + £~ (£ = e,fi). 

For the signal process, e + e~ — > pn, the Monte Carlo events are generated with an angular 
distribution of sin 2 #i(l + cos 2 9 V + sin 2 9 V cos(20i)), where 9 V is the angle between the p and 
the positron direction, and 9\ and 0i are the polar and azimuthal angles of the pion in 
the p helicity frame. A Monte Carlo sample of ip" ~^ vr + 7r _ 7r° is also generated with the 
same angular distribution, taking into account the available phase-space of the two-pion 
system. The generators include the effects of initial state radiation (ISR), and the pn or 
7T + 7r~7r° form factor varies as a function of s (s denotes the square of the center-of-mass 
energy), where a 1/s dependence is assumed. The generator PPCON is also used to study 
the interference between the continuum and resonant ip" amplitudes. In this generator, 
the branching fraction of ip" — > pit can be set to be any number between and 1, the 
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relative phase between the ip" strong and electromagnetic decay amplitudes can be set to 
be any possible value from —180° to 180° and the measurement of a Born (e + e~ — > pit) at 
y/s = 3.67 GeV measured by CLEO-c [4] is used to normalize the contribution of the 
continuum cross section at the ip" peak. 

Monte Carlo samples of Bhabha, dimuon, DD and inclusive hadronic events generated 
with Lundcrm [13] are used for the background study. 

II. EVENT SELECTION 

The final states of the study include two charged pions and one neutral pion which is 
reconstructed from two photons. Event selection includes photon identification and charged 
particle identification. 

A neutral cluster is considered to be a photon candidate when the deposited energy in the 
BSC is greater than 80 MeV, the angle between the nearest charged track and the cluster is 
greater than 16°, and the angle between the two nearest photons is larger than 7°. The first 
hit of the cluster is in the beginning 6 radiation lengths, and the angle between the cluster 
development direction in the BSC and the photon emission direction must be less than 37°. 
The number of photon candidates after the above selection is required to be two. 

For each charged track, the TOF and dE/dx measurements are used to calculate x 2 
values and the corresponding confidence levels to the hypotheses that the particle is a pion, 
a kaon or a proton (Prob w , Probx, Prob p ). A track is considered to be a pion when the 
confidence level of the pion hypothesis is greater than the confidence levels of the kaon and 
proton hypotheses. At least one charged track is required to be identified as a pion. 

For the decay channel of interest, the candidate events must satisfy the following selection 
criteria: 

1. An event is required to have only two oppositely charged tracks in the MDC, each with 
a good helix fit. The closest approach of the track to the interaction point is required 
to be within 2 cm in the transverse plane and within 20 cm in the beam direction, and 
the transverse momentum P xy > 0.06 GeV/c is used to remove the beam associated 
background. 

2. A four-constraint kinematic fit is performed under the hypothesis e + e~ — > jjn + n^. 
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and the confidence level of the fit is required to be greater than 1%. A four-constraint 
kinematic fit is also performed under the hypothesis of e + e~ — > r yyK + K~~, and x^-ymr < 
Xjjkk * s required to remove the K + K~n° events. 

3. To remove the di-muon background and the backgrounds produced by the ISR process 
e + e~ — > 7?//, with ip' — > neutral + J/ip, J ftp — > fi + {i~, two tracks should be in the 
| cos#| < 0.80 region (0 is the polar angle of the track in MDC) and at least one track 
is required to be in the coverage of muon counter, in which N h + + N^J < 3 is required. 
Here, N hlt is the number of muon counter layers with matched hits and ranges from 
to 3, indicating not a muon (0), a weakly (1), medianly (2), or strongly (3) identified 
muon track [14] . 

4. After the four-constraint kinematic fit, the energy of the higher momentum photon 
candidate is required to be less than 1.5 GeV to remove the p°(770) background 
produced by ISR. 

5. After the above selection, the radiative Bhabha background can still be seen clearly 
from Figure 1, where the dE/dx separation from pion hypothesis (XdE/dx) an< ^ the 
energy deposited in the BSC of the charged track (E B sc) are shown. The cluster of 
events in the top right corner are electron tracks and can be removed by requiring 
XdE/dx < —2Ebsc + 3 with high efficiency for the signal events. 

After applying all of the above selection criteria, the invariant mass distributions of the 
two photons after the kinematic fit are shown in Figures 2a and 2b for the ip" and the 
continuum data samples, respectively. It can be seen that there are 7r° signals in the ip" and 
y/s = 3.65 GeV data samples. 

Analyses using the Monte Carlo samples of Bhabha, dimuon, DD, inclusive hadronic 
events and ISR production of J / ip and ip' decays show that the background contaminations 
to the two-photon invariant mass spectra are small or in a random distribution. 

III. FITS TO THE INVARIANT MASS SPECTRA 

The invariant mass spectra of the two photons are fitted with the Monte Carlo simulated 
7T° invariant mass distribution (where the 7r° mass is fixed to the PDG value [15] in the 
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FIG. 1: Scatter plot of the dE/dx separation from the pion hypothesis versus the deposited energy 
in BSC of the charged track. The cluster at the top right corner is radiative Bhabha electrons. 
The events above the straight line are removed as Bhabha candidates. 
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FIG. 2: Invariant mass distributions of the two photons after final selection for (a) ip" data sample, 
and for (b) y/s = 3.65 GeV data. The histograms are data, and the curves show the best fits. 

simulation) for the signal and a 2nd-order polynomial for the background. The fits yield 
11.4 ± 4.7 and 10.0 ± 3.8 n°s for the ip" and the \fs = 3.65 GeV data samples respectively, 
and the corresponding signal significance are 3.1a and 4.6a. The fit results are shown in 
Figures 2a and 2b. The efficiency of detecting e + e~ — > it + it~it at ip" is (7.65 ± 0.12)% and 
that at y/s = 3.65 GeV is (7.88 ±0.12)% according to the Monte Carlo simulation by assum- 
ing a phase space distribution, where the errors are statistical due to limited statistics of the 
Monte Carlo samples. Here in the generator, only the continuum amplitude is considered. 

In order to evaluate the contributions of pn in the data samples, the Dalitz plots of the 
7T + 7r~7r° system are shown in Figures 3a and 3b after requiring the invariant mass of the two 



photons lies between 0.10 and 0.17 GeV/c 2 (about 2a around the n° nominal mass). The 
Monte Carlo predicted pn events are also shown in the plots as small black dots. No clear 
p(770) signal can be seen either in the ip" or in the continuum data sample. 
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FIG. 3: Dalitz plots of e + e~ — ► n + n~n° at (a) ip", and (b) y/s = 3.65 GeV after the final selection. 
The big black dots are for data, and the small black dots are the Monte Carlo simulated pn events. 
There are still non-7r° backgrounds in both plots for data. 



Since the pn signals are not significant, we try to set the upper limits for the e + e~ — > pn 
cross sections at ip" an d at y/s = 3.65 GeV. This is done by fitting the two photon invariant 
mass distributions (shown in Figure 4) after requiring the n + n~ or r yyn ± invariant mass in 
the range of 0.626 to 0.926 GeV/c 2 . The fits yield 2.9 ± 2.8 and 4.9 ± 2.7 n°s at f and 
y/s = 3.65 GeV respectively. The upper limits at 90% C. L. on the numbers of events are 5.1 
and 8.3 at ip" and y/s = 3.65 GeV respectively. Here the systematic errors are considered, 
see below. The number of pn events are overestimated since there are, in general, non-p7r 
contributions to events containing 7r°'s. The detection efficiency is (4.87 ± 0.09)% at ip" 
and (5.14 ± 0.09)% at y/s = 3.65 GeV from the Monte Carlo simulations assuming a pure 
continuum contribution. 



IV. SYSTEMATIC ERRORS 

Many sources of systematic error are considered in the cross section measurement. Sys- 
tematic errors associated with the efficiency are determined by comparing J/ip and ip' 
data and Monte Carlo simulation for very clean decay channels, such as J/ip — > n + n~n°, 
ip' — > n + n~J/ip, which allow the determination of systematic errors associated with the 
MDC tracking, trigger, kinematic fitting, photon identification efficiency, requirement of 
photon number and particle identification [16, 17]. The uncertainty due to the generator for 
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FIG. 4: Invariant mass distributions of the two photons of e + e~ — > pit at (a) ip" and (b) yfs = 
3.65 GeV after making the p mass requirement. The histograms are data, and the curves show the 
best fits. 

the it + it~it° mode is estimated by the efficiency difference between the phase space and the 
pit generators. The systematic errors from the background estimation are determined by 
comparing the fitting result between the 2nd-order polynomial background and a lst-order 
polynomial background and different fitting ranges in the fits of the two-photon invariant 
mass spectra. Uncertainty of the integral luminosity of the data sample is also a source of 
the systematic error. 

All the sources considered are listed in Table I. The total systematic errors for 



pit are 20.8% and 16.4% for the ip" and the continuum data, respectively, and 



for e + e — > it + it it the corresponding numbers are 23.8% and 19.3% for ip" and continuum 
data, respectively. 



V. RESULTS AND DISCUSSION 

We give the measurement of the cross section of e + e~ — > it + it~it° and the upper limit of 
the cross section of e + e~ — > pit at ip" and y/s = 3.65 GeV. For the process e + e~ — > it + it~it°, 
the Born order cross section is calculated with 

a fe e — > it it it = — — — — , 

and for e + e~ — > pit, the upper limit of the cross section is calculated with 

jvrobs 



e ■ £ ■ B(it° ^ + 5) 
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TABLE I: Systematic errors for pir and tt + tt tt° cross section measurements at tp" and y/s = 
3.65 GeV (%). 



Source 


y/s = 3.65 GeV 


y/' 


Monte Carlo statistics 


1.7 


1.7 


Trigger 


1.0 


MDC tracking 


4.0 


Kinematic fit 


6.0 


Photon efficiency 


4.0 


Number of photons 


2.0 


Particle ID 


negligible 


Background 


13.4 


18.7 


Luminosity 


3.7 


3.0 


Generator (for 7r + 7r~7r° only) 


10.2 


11.5 


Total [pir (vr+vr-vr )] 


16.4 (19.3) 


20.8 (23.8) 



Here the initial state radiative correction factor, 1 + 5, and the efficiency, e, are ob- 
tained from the Monte Carlo simulation assuming a pure continuum contribution. Using 
the numbers obtained from the above analysis (listed in Table II), one gets 

a B (e + e- -> tt+tt^tt ) = (8.4 ± 3.5 ± 2.0) pb, 
a B (e + e~ -> pir) < 6.0 pb 

at ip" and 

a(e + e- -> tt+tt-tt ) = (19.3 ± 7.3 ± 3.7) p6, 
cr(e + e _ — > p7r) < 25 p6 

at v^s = 3.65 GeV, where the first errors are statistical and the second ones are systematic, 
and the upper limits are at 90% C. L. 

The upper limit of e + e _ — > pn cross section at y/s = 3.65 GeV is consistent with the 
measurement at y/s = 3.67 GeV by CLEO-c [4] and the calculation from Ref. [7]. The upper 
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TABLE II: Numbers used in the calculations of e + e — > tt + tt tt and /97r cross sections. 



Decay Channel 


v/i = 3.65 GeV 




7T+7T 7T° 


pir 


7T+7T 7T° 


pir 


N obs {N^ b £) 


10.0 ±3.8 


4.9 ±2.7 (8.3) 


11.4 ±4.7 


2.9 ± 2.8 (5.1) 


e(%) 


7.88 ±0.12 


5.14 ±0.09 


7.65 ±0.12 


4.87 ±0.09 


1 + S 


1.033 


1.025 


1.033 


1.026 


Lipb- 1 ) 


6.42 ±0.24 


17.3 ±0.5 


B(ir° -» 77) [15] 


0.988 



limit of o"(e + e~ — > p7r) at ip" is a little lower than the upper limit of 6.3 pb determined by 
Mark III [8] and is lower than the measurement of <r(e + e - — > pn) = 8.0l];' 4 ± 0.9 pb off the 
ip" peak at y/s = 3.67 GeV by CLEO-c [4], which indicates that there must be a non-zero 
■ip" — > pn amplitude at the ip" energy. 

The measurement of <x(e + e~ — > pit) at -0" also supports the postulation in Ref. [7] that 
the relative phase between the strong and electromagnetic decays of tp" into light hadrons is 
around —90°. In this scheme [7], the number of observed pn events at the ip" peak depends 
on both the ifj" — > pn branching fraction B(ip" — > pir) and the relative phase, 0, between the 
ip" strong and electromagnetic decay amplitudes. Using the measurement of pn at the ip" 
peak in this experiment, and the e + e~ — > pit cross section at the continuum by the CLEO-c 
experiment [4], a 2-dimensional scan indicates that the physical region of B(ip" — > pn) and 
4> is restricted in the hatched area as shown in Figure 5 at 90% C. L. If the correlation 
between B(ip" — > pit) and is neglected, one gets B(ip" — > pn) e (6.0 x 10" 6 , 2.4 x 10~ 3 ), 
and G (—150°, — 20°) at 90% C. L. This result is consistent with the calculation in 
Refs. [6, 7] that the ip" — > pn branching fraction is at the 10~ 4 level, and supports the 
explanation of the " pn puzzle" observed between J/tp and ip' decays by the S- and D-wave 
mixing model. We also expect CLEO-c, with a few fb~ l of ip" data [18], and BESIII, with 
even more ip" data [19], to be able to produce tighter constraints on B(ip" — > pir) and 0, 
and give a better test of this scenario. 
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-ISO -160 -140 -120 -100 -80 -60 -40 -20 
Phase between \\f" Strong and EM decays 

FIG. 5: Restriction on B(ip" — > pir) and <p from the measurement of pir at the ip" peak in this 
experiment. The hatched area indicate the physical region at 90% C. L. 

VI. SUMMARY 

The processes e + e~ — > it + it^it and pn are searched for at ip" and y/s = 3.65 GeV. 
We observe e + e~ — > 7r + 7r~7r° signals at the 3.1<r and 4.6<7 levels for y/s = 3.773 GeV and 
y/s = 3.65 GeV, respectively. No significant pir signal is observed, and the upper limit of 
the e + e" — > pn cross section at ip" is measured to be 6.0 pb at 90% C. L. assuming no 
ip" contribution. Considering the interference between the continuum amplitude and the 
ip" resonance amplitude, the branching fraction of ip" decays to pir is determined to be 
B(ip" -> pvr) G (6.0 x 10~ 6 , 2.4 x 10~ 3 ) at 90% C. L. 
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